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ABSTRACT 

X-ray observations of galaxy clusters potentially provide powerful cosmological probes if systematics due 
to our incomplete knowledge of the intracluster medium (ICM) physics are understood and controlled. In this 
paper, we present mock Chandra analyses of cosmological cluster simulations and assess X-ray measurements 
of galaxy cluster properties using a model and procedure essentially identical to that used in real data analysis. 
We show that reconstruction of three-dimensional ICM density and temperature profiles is excellent for relaxed 
clusters, but still reasonably accurate for unrelaxed systems. The total ICM mass is measured quite accurately 
(<6%) in all clusters, while the hydrostatic estimate of the gravitationally bound mass is biased low by about 
5%-20% through the virial region, primarily due to additional pressure support provided by subsonic bulk 
motions in the ICM, ubiquitous in our simulations even in relaxed systems. Gas fraction determinations are 
therefore biased high; the bias increases toward cluster outskirts and depends sensitively on its dynamical 
state, but we do not observe significant trends of the bias with cluster mass or redshift. We also find that 
different average ICM temperatures, such as the X-ray spectroscopic r spec and gas-mass-weighted T m$ , are 
related to each other by a constant factor with a relatively small object-to-object scatter and no systematic trend 
with mass, redshift or the dynamical state of clusters. We briefly discuss direct applications of our results for 
different cluster-based cosmological tests. 

Subject headings: cosmology: theory-clusters: formation- methods: numerical 



1. INTRODUCTION 

X-ray observations of clusters of galaxies can poten- 
tially provide plethora of useful cosmological information. 
The cluster-based cosmological tests include cluster number 
counts based on the temperature function, baryon fraction and 
X-ray and SZ Hubble constant measurements. These cosmo- 
logical tests provide potentially very powerful constraints on 
the matter density, £Im, Qa, erg, and the equation of state of 
dark energy, w. In all of these cosmological tests, the key 
observational ingredients are the gas mass (M gas ), total mass 
(M to t), gas mass fractions (/ gas =M gas /M tot ), and average clus- 
ter temperature ((T)). In the era of precision cosmology, it is 
paramount to achieve accurate measurements of the key clus- 
ter parameters, check how they depend on various simplifying 
assumptions, and control systematics due to our incomplete 
knowledge of the ICM physics. 

There are many important simplifying assumptions used in 
deriving cluster properties from X-ray data. For example, it is 
usually assumed that the ICM density is a function of radius 
only and does not have small-scale substructure. However, 
substructure in clusters is ubiquitous and this biases the X-ray 
derived M gas high (Mathiesen et al. 1999). Also, the small- 
scale clumps are often associated with subhalos and thus have 
lower temperature. Neglecting this can bias the average spec- 
troscopic temperature of the cluster as a whole (Mathiesen 
& Evrard 2001; Mazzotta et al. 2004) which could poten- 
tially lead to biases in constraints on cosmological parame- 
ters (Rasia et al. 2005). Cluster total masses are often esti- 
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mated assuming the hydrostatic equilibrium for the ICM. This 
may underestimate the true gravitational bound mass if non- 
thermal pressure support is present. For example, turbulent 
gas motions can provide about 10%-20% of the total pres- 
sure support even in relaxed clusters and hence bias the hy- 
drostatic estimate (Schuecker et al. 2004; Faltenbacher et al. 
2005; Rasia et al. 2006; Lau et al. 2006). Spherical symmetry 
is another commonly used assumption useful in solving for 
three-dimensional (3D) physical cluster properties from two- 
dimensional (2D) observed quantities, but could bias results 
if clusters are aspherical. 

X-ray observations with Chandra and XMM-Newton en- 
able us to study properties of the ICM with unprecedented 
detail and accuracy and provide important handles on the ICM 
modeling and associated systematics. Their superb spatial 
resolution and sensitivity enable accurate X-ray brightness 
and temperature measurements at a large fraction of the clus- 
ter virial radius and also make it simple to detect most of the 
small-scale X-ray clumps. Despite this recent observational 
progress, the biases in the determination of the key cluster 
properties remain relatively uncertain. The main obstacle is 
that because of their unrivaled statistical accuracy, the X-ray 
results cannot be contrasted against any other independent ob- 
servational techniques. In this study, we attempt to check the 
validity of the X-ray analyses, with a specific focus on the 
analyses of high-resolution Chandra data, by using mock ob- 
servations of clusters derived from cosmological simulations, 
for which the true answers are known. 

Our simulations properly treat both collisionless dynamics 
of dark matter and stars and gasdynamics and capture a va- 
riety of physical phenomena from the nonlinear collapse and 
merging of dark matter to shock-heating and radiative cool- 
ing of gas, star formation, chemical enrichment of the ICM by 
supernova and energy feedback. These simulations can there- 
fore be used to test observational biases due to incomplete re- 
laxation of gas, the dynamical state of a cluster, substructure, 
or non-isothermality. As we show in the forthcoming papers, 
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the current simulations match the observed ICM profiles out- 
side cluster cores (Nagai et al. 2006) and the global gas mass 
fraction (Kravtsov et al. 2006a) and are therefore sufficiently 
realistic for a purpose of the current study. 

The mock observations of the simulated clusters are gen- 
erated by exactly reproducing the response properties of the 
Chandra telescope. The analysis of the mock data is identi- 
cal to that used for real Chandra observations by Vikhlinin 
et al. (2005, 2006). The comparison of the true and derived 
cluster properties provides an assessment of any biases intro- 
duced by the X-ray analysis. Some of these biases are poten- 
tially redshift-dependent (e.g, those related to higher merger 
rate or decreased ability to detect small-scale substructures in 
the X-ray data for high-z systems). To check for any redshift 
dependence in such biases we use the simulation outputs at 
z = and 0.6. Our results indicate that the X-ray analysis 
provides accurate reconstruction of the 3D properties of the 
ICM. The strongest biases we find are those in the hydrostatic 
mass estimates. They are related to physics explicitly miss- 
ing from the hydrostatic method (e.g., turbulence), and not to 
deficiencies of the X-ray analysis. 

2. COSMOLOGICAL CLUSTER SIMULATIONS 

In this study, we analyze high-resolution cosmological 
simulations of 16 cluster-sized systems in the flat ACDM 
model: Q m = 1 - D. A = 0.3, D. b = 0.04286, h = 0.7 
and erg = 0.9, where the Hubble constant is defined as 
100/z km s _1 Mpc~', and an <x 8 is the power spectrum normal- 
ization on an 8/z~' Mpc scale. The simulations were done with 
the Adaptive Refinement Tree (ART) A^-body+gasdynamics 
code (Kravtsov 1999; Kravtsov et al. 2002), a Eulerian code 
that uses adaptive refinement in space and time, and (non- 
adaptive) refinement in mass (Klypin et al. 2001) to reach 
the high dynamic range required to resolve cores of halos 
formed in self-consistent cosmological simulations. The sim- 
ulations presented here are discussed in detail in Kravtsov 
et al. (2006b) and Nagai et al. (2006) and we refer the reader 
to these papers for more details. Here we summarize the main 
parameters of the simulations. 

High-resolution simulations were run using a uniform 128 3 
grid and 8 levels of mesh refinement in the computational 
boxes of 120 h~ x Mpc for CL101-107 and 80 h~ x Mpc for 
CL3-24 (see Table 1). These simulations achieve a dynamic 
range of 32768 and a formal peak resolution of » 3.66 h~ x kpc 
and 2.44 hr x kpc, corresponding to the actual resolution of 
KlhT x kpc and 5 h~ x kpc for the 120 and 80 h~ x Mpc boxes, 
respectively. Only the region of ~ 3 - 10 /T 1 Mpc around the 
cluster was adaptively refined, the rest of the volume was fol- 
lowed on the uniform 128 3 grid. The mass resolution, m pa rt, 
corresponds to the effective 512 3 particles in the entire box, 
or the Nyquist wavelength of A Ny = 0.469 h~ x and 0.3 12 ft -1 
comoving Mpc for CLIO 1-1 07 and CL3-24, respectively, or 
0.018 h~ x and 0.006 hr x Mpc in the physical units at the initial 
redshift of the simulations. The dark matter particle mass in 
the region around the cluster was 9.1 x 10 8 h~ x M© for CL101- 
107 and 2.7 x 10 8 /i 1 M for CL3-24, while other regions 
were simulated with lower mass resolution. 

The Af-body+gasdynamics cluster simulations used in this 
analysis include collisionless dynamics of dark matter and 
stars, gasdynamics and several physical processes critical to 
various aspects of galaxy formation: star formation, metal 
enrichment and thermal feedback due to supernovae Type II 
and Type la, self-consistent advection of metals, metallicity 
dependent radiative cooling and UV heating due to cosmo- 



TABLE 1 

Simulated cluster sample of the CSF run at z=0 



Name 


(A - 'Mpc) 


M gas 

5QQc 

(/r'l0 fi M o ) 


ilftot 
5QQc 

(/r'lO^Mo) 


(keV) 


CL101 


1.160 


8.19 


9.08 


8.7 


CL102 


0.978 


4.83 


5.45 


5.6 


CL103 


0.993 


4.93 


5.71 


4.7 


CL104 


0.976 


5.17 


5.39 


7.7 


CL105 


0.943 


4.73 


4.86 


6.2 


CL106 


0.842 


3.18 


3.47 


4.3 


CL107 


0.762 


2.17 


2.57 


4.0 


CL3 


0.711 


1.92 


2.09 


3.7 


CL5 


0.609 


1.07 


1.31 


1.9 


CL6 


0.661 


1.38 


1.68 


2.0 


CL7 


0.624 


1.22 


1.41 


1.9 


CL9 


0.522 


0.74 


0.82 


1.1 


CLIO 


0.487 


0.43 


0.67 


2.4 


CL11 


0.537 


0.78 


0.90 


3.4 


CL14 


0.509 


0.62 


0.77 


3.0 


CL24 


0.391 


0.26 


0.35 


1.5 



"Spectral temperature measured in the shell of [0.15,l]r5oo c - 



logical ionizing background (Haardt & Madau 1996). The 
cooling and heating rates take into account Compton heating 
and cooling of plasma, UV heating, and atomic and molec- 
ular cooling, and are tabulated for the temperature range 
10 2 < T < 10 9 K and a grid of metallicities, and UV inten- 
sities using the Cloudy code (ver. 96b4; Ferland et al. 1998). 
The Cloudy cooling and heating rates take into account metal- 
licity of the gas, which is calculated self-consistently in the 
simulation, so that the local cooling rates depend on the local 
metallicity of the gas. Star formation in these simulations was 
done using the observationally-motivated recipe (e.g., Kenni- 
cutt 1998): p. = Pg a 5 s /f., with U = 4 X 10 9 yrs. The code 
also accounts for the stellar feedback on the surrounding gas, 
including injection of energy and heavy elements (metals) via 
stellar winds, supernovae, and secular mass loss. 

These simulations therefore follow the formation of galaxy 
clusters starting from the well-defined cosmological initial 
conditions and capture the dynamics and properties of the 
ICM in a realistic cosmological context. However, some po- 
tentially relevant physical processes, such as AGN bubbles, 
magnetic field, and cosmic rays, are not included. Therefore, 
the simulated ICM properties are probably not fully realis- 
tic in the innermost cluster regions. Moreover, the gas in the 
simulations is treated as an ideal inviscid fluid with a small 
amount of numerical viscosity, and it remains unclear to what 
extent the level of ICM turbulence found in the simulations 
and discussed below applies to real clusters. Despite these 
limitations, the current simulations reproduce the observed 
ICM profiles outside cluster cores (Nagai et al. 2006) and the 
global gas mass fraction (Kravtsov et al. 2006a) and are there- 
fore sufficiently realistic for a purpose of the current study. 

Our simulated sample includes 16 clusters at z = and 
their most massive progenitors at z = 0.6. The properties of 
simulated clusters at z = are given in Table 1 . The masses 
are reported at the radius rsm c enclosing overdensities with re- 
spect to the critical density at the redshift of the output (below, 
we also use a higher overdensity level, 2500). 

The mass measurement biases obviously need to be consid- 
ered separately for dynamically relaxed and non-relaxed clus- 
ters. Our relaxed subsample is identified based on the overall 
structural morphology of their Chandra images, mimicking 
the procedure used by observers. Specifically, we visually 
examine mock 100 ksec images and identify "relaxed" clus- 
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Fig. 1. — Mock Chandra images of the z = simulated clusters. The left and right panels show the images of the relaxed cluster CL104 and unrelaxed CL101, 
respectively. The detectable extended X-ray sources (other than the cluster itself) are detected and masked out from the analysis (red ellipses). The physical size 
of the images is 5 Mpc. 



ters as those with regular X-ray morphology and no secondary 
maxima and minimal deviations from elliptical symmetry. By 
contrast, "unrelaxed" clusters are those with secondary max- 
ima, filamentary X-ray structures, or significant isophotal cen- 
troid shifts. The typical examples of systems classified as re- 
laxed or unrelaxed are shown in Figure 1 . 

3. MOCK CHANDRA ANALYSES 

3.1. Simulating the Mock Chandra Data 

We create X-ray flux maps of the simulated clusters viewed 
along three orthogonal projections. The flux map is com- 
puted by projecting X-ray emission arising from hydrody- 
namic cells enclosed within 3 R vn of a cluster along a given 
line-of-sight. We compute the X-ray spectrum arising from 
the /-th hydrodynamical cell with a cell volume V, as 



h,i = n ei n pi l^ E {JuZuz)Vu 



(1) 



where n e j, n p j, T,, and Z, are the electron and proton densities, 
gas temperature, and metal abundance of the hydrodynamic 
cell. We compute the X-ray plasma emissivity, A£(r,-,Zi, z), 
using the mekal code (Mewe et al. 1985; Kaastra & Mewe 
1993; Liedahl et al. 1995) with the relative heavy element 
abundance from Anders & Grevesse (1989). We do not com- 
pute emission from the gas with T < 10 6 K (0.086 keV) be- 
cause it is below the lower limit of the Chandra bandpass. 
The plasma spectrum is multiplied by the photoelectric ab- 
sorption corresponding to the hydrogen column density of 
nil = 2xl0 20 cirT 2 . The resulting map consists of 1024x1024 
2D image at 218 energy bins ranging from 0.1 to lOkeV. The 
pixel size is 4.88 h' 1 kpc for CL101-CL107 and 2.44 ft" 1 kpc 
for CL3-CL24. The map size is therefore 5h ~' Mpc and 
2.5 hT l Mpc, respectively. Throughout this paper, we assume 
the cluster redshift of z Q b s = 0.06 for the 2=0 sample and 
^obs = 0.6 for the z = 0.6 sample. 

Next, we simulate mock Chandra data. We convolve the 
emission spectrum with the response of the Chandra front- 
illuminated CCDs and draw a number of photons at each po- 



sition and spectral channel from the Poisson distribution. We 
simulate two sets of mock Chandra data. The first set has 
an exposure time of 100 ksec (typical for deep observations) 
and includes a background with the intensity corresponding 
to the quiescent background level in the ACIS-I observations 
(Markevitch et al. 2003). From these data, we generate im- 
ages in the 0.7-2 keV band and use them to identify and 
mask out from the further analysis all detectable small-scale 
clumps, as routinely done by observers. Our clump detection 
is fully automatic and based on the wavelet decomposition al- 
gorithm described in Vikhlinin et al. (1998). Detection thresh- 
olds were chosen to allow 3-4- false detections in each image. 
The typical limiting flux for detection of compact extended 
sources is ~ 3 x 10~ 15 erg s ctrT 2 in the 0.5-2 keV band; 
this corresponds to a luminosity of ~ 1.5 x 10 42 erg s^ 1 for 
clusters at z D b s = 0.06. Detected clumps are indicated by red 
ellipses in the images presented in Fig.l. Note that the M gas 
associated with the excluded clumps is no more than a few 
percent of the total enclosed M gas within rsooc at z — 0.06. 
At z = 0.6, considerably smaller numbers of clumps are de- 
tected 5 , and the clump removal has very little effect on the 
analysis of high-z clusters. 

Once the clump detection is done on the first set, all further 
analysis is performed using a second set of photon maps gen- 
erated for very long exposures yielding a 10 6 - 10 7 photons 
outside a cluster core region in each of the simulated clusters. 
The exposures are artificially long by design as we are inter- 
ested in exploring intrinsic limitations of the X-ray analysis, 
not the statistical errors due to Poisson noise in a particular 
choice of short exposure. Also, we ignore further compli- 
cations present in reduction of real Chandra data, including 
background subtraction and spatial variations of the effective 
area (i.e., we assume that accurate corrections for these effects 
can be applied to the real data and any associated uncertainties 

5 Only a few clumps are detected in massive clusters with r spcc > 3 keV 
and none in less massive systems. 
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are included in the reported measurement errors). 

3.2. Analysis of Mock Chandra Data 

Using the simulated mock data as an input, we repeat 
the analysis applied by Vikhlinin et al. (2006) and Kotov & 
Vikhlinin (2006) to real Chandra observations of relaxed clus- 
ters at z ~ and 0.5. The purpose of this analysis is to re- 
construct spherically averaged ICM density and temperature 
profiles, and to estimate total cluster mass assuming hydro- 
static equilibrium. The entire procedure and justification of 
the choices made is extensively discussed in Vikhlinin et al. 
(2006). Below, we briefly outline the essential steps. 

The main observational input is azimuthally averaged X-ray 
surface brightness and projected gas temperature and metal- 
licity profiles. The temperature and metallicity profiles are 
obtained from a single-temperature fit to the spectra in the 
0.5-10 keV extracted in annuli with r out /r m = 1.25 centered 
on a main cluster X-ray peak. 

The surface brightness profiles are extracted in the 0.7- 
2 keV band in narrow concentric annuli (r out /r m = 1.1). Us- 
ing the effective area as a function of energy and observed 
projected temperature and metallicity at each radius, we then 
convert the observed Chandra count rate in the 0.7-2 keV 
band into the projected emission measure, EMM = Jn e n p dl. 

Three-dimensional profiles of p gds (r) and r gas (r) are recon- 
structed by using analytic 3D models with great functional 
freedom, projecting them along a line of sight, and fitting the 
projected models to observed profiles. The gas density model 
we use is given by 



(r/r c r a 



1 



(i + r 2ir 2 c fP- a i 2 (i + ry/r s yy'y' 



(2) 



This expression is a modification of the j6-model (Cavaliere 
& Fusco-Femiano 1978) that allows independent modeling of 
changes in the gas density slopes in cluster cores, outskirts, 
and intermediate radii. Vikhlinin et al. (2006) used an addi- 
tional /3-model component with small core-radius to further 
increase functional freedom in a very central region. We do 
not use this feature and instead focus on modeling the surface 
brightness in a region at r > 0.06 x rsooc- 

The analytic expression of a 3D temperature profile was de- 
signed to model general features of observed projected tem- 
perature profiles (Vikhlinin et al. 2005). The model consists 
of two terms: 

T 3D (r) = t(r) x f cool (r). (3) 

The first term, t(r), describes a temperature decline at large 
radii by a broken power law profile with a transition region, 



t(r) = T 



(r/r t y a 



(l+(r/r t ) b y/f 



(4) 



The second term, t CO oi(r), is designed to model the temperature 
decline in the central region affected by radiative cooling 



fcooi(r) = (x + T min /T )/(x +1), x = (r/r cooi ) aa 



(5) 



To fit the observed projected temperature profile, the model 
is projected along a line of sight. This projection requires 
proper weighting of multiple temperature components, and 
we use the algorithm that very accurately predicts a single- 
temperature fit to multi-component spectra in a wide range of 
temperatures recently proposed by Vikhlinin (2006). The ob- 
served projected temperature profiles are fit in the radial range 
between O.lrsooc and r out = max(1.5r5oo c , 2/3r max ), where 
'"max is a half of the mock Chandra image size. The inner 



radial range is set to exclude central components with unreal- 
istic X-ray properties found in most of our simulated clusters. 
The exact choice of r m [ n is unimportant because we are pri- 
marily interested in the cluster properties at large radii. The 
choice of the outer radius is motivated by a desire to ensure 
accurate measurements of the temperature gradient and hence 
the M tot at rsooc an d that the fits in the cluster outskirts are not 
affected by projection effects. 

Several examples of the observed surface brightness and 
projected temperature profiles are shown in Figures 2 and 3 
along with their best fit models. We find that the models of 
eq.(2) and (3) usually provide an excellent description of the 
data, except for very irregular clusters such as that in Figure 3 
(its X-ray image is shown in the right panel of Figure 1). 

Given 3D models for the gas density and temperature pro- 
files, the hydrostatic estimate of total cluster mass within a 
radius r is given by 



M(r) = — 



T(r)r Id log p g d log T 



(6) 



fi m p G \ d log r d log r ) ' 

where fi = 0.592 is mean molecular weight of the fully ion- 
ized H-He plasma used in the simulations. Given M(r), we 
compute the total mass at several overdensity levels, A, rel- 
ative to the critical density at the cluster redshift, by solving 
equation 

M A (r A ) = A4/3 nr{p c (z). (7) 

The ICM mass profile is obtained directly from the best-fit 
number density profile (eq.2), p g = 1.236 m p (n p n e )^ 2 . 

Finally, we compute average temperatures for each clus- 
ter using different weightings of the best-fit 3D temperature 
models and observed projected profiles. In this study, we con- 
sider several definitions of average temperatures. The first is 
a X-ray spectral temperature, r spec , a value derived from a 
single-temperature fit to the integrated cluster spectrum ex- 
cluding the core and detectable small-scale clumps. We also 
compute a gas-mass-weighted temperature, T mg , obtained by 
weighting of the reconstructed temperature profile with the 
derived gas density profile. We also consider an emission- 
weighted temperature, T^ 1 , computed by weighting T(r) 
with n p n e T 1 ^ 2 ; it is not directly observed but often used in 
the theoretical work. Similarly, we compute a true gas-mass- 
weighted temperature T^ 1 from the simulations. The super- 
script "SIM" indicates that the quantity is computed directly 
from the 3D gas properties in the simulations without exclud- 
ing the Chandra detectable clumps, as opposed to from the 
Chandra analysis fits, throughout this work. All these average 
temperatures are computed excluding the central 0. 15 rsooc be- 
cause temperatures at these radii can be strongly affected by 
radiative cooling and thus not directly related to the depth 
of the cluster potential well. Also, because of the complex 
physics in the innermost cluster regions the simulated ICM 
properties are probably not fully realistic. The outer radius of 
integration is set to rsooc, which approximately corresponds 
to the boundary within which the currently observed cluster 
temperature profiles are reliable (Vikhlinin et al. 2005). How- 
ever, we also consider a spectral temperature r spe c[0.5r5ooc] 
measured within a smaller aperture of 0.15 < r/rsooc < 0.5; 
this definition may be more practical than r spec integrated out 
to rsooc for observations of limited statistical quality or spa- 
tial coverage. Note that T svec is equivalent to T x defined in 
Kravtsov et al. (2006b). 

4. RESULTS 
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Fig. 2. — The mock Chandra analyses of a typical relaxed cluster, CL7, at z = 0. The cluster has r^ooc = 891kpc, M500C = 1 -41 X 10 14 /i _1 M Q and 
^sooc ~ 0.5r500c- The upper-left panel shows the observed emission measure profile (« surface brightness). The measured projected temperature profile is shown 
in the upper-middle panel. In these panels, the dashed lines indicate best-fit projected profiles plotted for the radial range where the fits were made. The lower-left 
and -middle panels show comparisons of the best-fit 3D model profiles and true gas density and temperature profiles. Note that our model recovers well both the 
projected profiles and the actual 3D profiles. In the upper-right panel, we compare the derived and true M tot (upper lines) and M gas profiles (lower lines). The 
reconstructed M gas profile dot-dashed line) is accurate to a few percent in the entire radial range shown. The hydrostatic M tot estimate (dotted line), on the other 
hand, is biased low by about 5%-10% in the radial range, [0.2, l.Ojrsooc. The lower-right shows that measured cumulative / eas is biased high by ~ 10% in the 
radial range of [0.2, 1.0]r5oo c for this cluster, and it is primarily due to the bias in the hydrostatic mass estimate. 



4. 1 . Profiles of Relaxed clusters 

Figure 2 shows the profile analysis of one of the relaxed 
clusters, CL7, at z = 0. The upper-left panel shows the 
emission measure (EMM) profile constructed from the mock 
Chandra photon map. The dashed line is a best-fit model in 
the radial range between 0.06r5Qo c and 2r 5 oo c - In this radial 
range, the EMM profile is very smooth and the model pro- 
vides an excellent fit. We do not model an additional compo- 
nent associated with dense and cool gas in the inner ss 50 kpc. 
Note, however, that the M gas (<50 kpc) is only a small fraction 

2%) of the M gas (< r 50 o c ). 

The upper-middle panel shows the projected temperature 
profile obtained from the X-ray spectral fitting at each pro- 
jected annulus. Data points are plotted as a function of X- 
ray counts weighted radius, and the error bars in the hor- 
izontal axis indicate a radial bin size within which the X- 
ray spectrum was extracted. The best-fit profile, indicated 
by the dashed line, is a very good description of the data in 
0.1 < r/r5oo c < 1.5, the radial range where the fit was done. 

In the lower-left and lower-middle panels, we compare the 
best-fit model profiles from the mock Chandra analysis to the 
true 3D gas density and temperature profiles measured in the 
simulation. For the true profiles, we compute a density profile 
of hot, X-ray emitting gas with T > 10 6 K and a gas-mass- 
weighted temperature. Note that a volume-averaged tempera- 



ture is quite similar to the gas-mass-weighted average shown 
here. For this cluster, both 3D gas density and temperature 
profiles are recovered very accurately. 

The upper-right panel shows the comparisons of the M tot 
and M gas profiles derived from the mock Chandra analysis 
and their respective true profiles. The comparisons show that 
the M gas profile is unbiased, and it is recovered to better than 
2% in 0.1 < r/rsooc < 2, nearly the entire radial range where 
the EMM profile was modeled. The M lot derived from hy- 
drostatic mass modeling, on the other hand, is biased low by 
about 5%-10% in the radial range of 0.2 < r/rsooc < 1- The 
figure also shows that the measured / gas is biased high, on av- 
erage, by about 10% at all radii, mainly due to the bias in the 
hydrostatic estimate of M tot . 

These results are representative of those for other clusters 
in the relaxed cluster samples at both z = and 0.6. The 3D 
gas density and temperature profiles are recovered remarkably 
well, and the M gas profiles are accurate to a few per cent at all 
radii. The M tot , on the other hand, is typically biased low by 
about 5%-20% throughout the clusters. One exception is the 
M tot measurement for one the most relaxed clusters, CL104, 
at z = 0, which is nearly unbiased. This system has not expe- 
rienced a major merger for several dynamical time 6 Gyr). 

4.2. Profiles of Unrelaxed Clusters 
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Fig. 3. — Same as Fig. 2, but for one of the uiirelaxed clusters, CL101 
at 2 = 0. The emission measure and projected temperature profiles exhibit 
bumps and wiggles. The derived M eas and M tot profiles are thus less accurate 
than those in relaxed systems. However, the fits to the projected data are still 
good, highlighting the flexibility of our gas density and temperature models. 



Accurate M tot estimates in unrelaxed clusters cannot be ob- 
tained from the hydrostatic equilibrium assumption. Recon- 
struction of 3D gas mass and temperature profiles from pro- 
jected X-ray data can be done only approximately in such sys- 
tems because of deviations from spherical symmetry. How- 
ever, many applications (e.g., cluster mass function estimates) 
require analyses of full statistical samples that include both re- 
laxed and unrelaxed systems. It is therefore important to ask 
how well the observational analysis procedures can recover 
properties of the ICM for unrelaxed systems. 

Unrelaxed clusters are a diverse population with a wide 
range of dynamical states. For the purpose of illustration, 
we show the profile analysis for a cluster undergoing a ma- 
jor merger, CL101 at z — (Fig. 3). This example highlights 
salient features of unrelaxed systems. Compared to relaxed 
clusters, their EMM and projected temperature profiles are 
less smooth and exhibit more pronounced bumps and wiggles. 
These features are not recovered by our smooth 3D model, 
even though the overall trends in the projected temperature 
and especially in the EMM are reproduced reasonably well. 
Recovery of the 3D gas density and temperature profiles is 
much less accurate than for relaxed systems. In the case of 
CL101, for example, the derived M gas profile is biased high 
by about 5%-10% at all radii, and the 3D temperature pro- 
file is biased low at r > O.lrsooc- If the hydrostatic equation 
is blindly applied to the best-fit models, the derived M tot is 
biased low by 25% at r^ooc- 

4.3. Gas Mass, Total Mass and Gas Mass Fractions 

In this section, we assess biases in the X-ray estimates of 
the enclosed M toU M gas , and / gas . The bias is defined as a 
fractional difference between estimated mass, M est , derived 
from the mock Chandra analysis and true 3D mass, M tme , 
measured directly in simulations, bias = (M est - M tme )/M tme . 
Figure 4 shows the results for all simulated clusters at z — 
and at radii enclosing two different overdensities: and 
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Fig. 4. — Bias in the M tot , Mga 8 , and / gas measured within r!™? (left) and 
r 2500c (rfshf) f° r simulated clusters at z = 0. The bias is defined as, e.g., 
(M:st - Mtrue)/Mtrue, where both M est and M trlle are at the same radius. Each 
cluster is viewed along three orthogonal projections, and clusters with relaxed 
and unrelaxed morphologies are indicated with filled and open symbols, re- 
spectively. 



r 500c' B° m estimated (M est ) and true masses (M tme ) are mea- 
sured in the same physical region enclosed within "true" radii 
('"true) measured directly from simulations. Note also that 
the Mg as is computed by excluding the mass enclosed within 
0.075r"g ( , c 6 (from both true and estimated masses), while the 
M tot is obtained without excluding the central region. The 
/gas = M g . d!i /M tot is therefore a fraction computed using these 
two quantities measured in a slightly different radial range. 
Table 2 also summarizes average biases and scatter in the 
mass estimates at both z — and 0.6. 

The upper panels in Figure 4 show that the hydrostatic mass 
estimate is typically biased low. It is underestimated on aver- 
age by 12% at rg , ( , 0c and 16% at r^ )c for all clusters at z = 0. 
The bias in the hydrostatic mass is smaller in the inner re- 
gion and it increases toward cluster outskirts. The bias is also 
smaller in relaxed systems than in unrelaxed systems; for ex- 
ample, biases are 8% and 15% at rS^, for the relaxed and 
unrelaxed samples, respectively. Note also that the scatter 
is very small inside ^socic °^ t ^ le re l axe d clusters. These re- 
sults indicate that the hydrostatic condition is best realized in 



The choice is made to ensure that the central, dense and cool gas compo- 
nent, associated with the ISM of the central cluster galaxy, is excluded from 
the M gas estimates of all clusters. This gas contributes < 2% of the total Mg as , 
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the inner region of relaxed clusters, while the deviation from 
the hydrostatic equilibrium become more prominent in clus- 
ter outskirts and/or unrelaxed systems, as expected. Similar 
results are obtained for clusters at high redshift. 

The middle panels in Figure 4 illustrate that the X-ray gas 
mass determinations are remarkably accurate and robust. The 
M gas measurements are accurate to better than a few percent 
for both relaxed and unrelaxed clusters and independent of 
redshift. For the relaxed cluster samples, the accuracy of the 
M gas measurements is as good as 1% at rSjjL with 3% 1 <x 
scatter. Although the effect is small (< 6%), the M gas is biased 
high in the outskirts of the unrelaxed systems. 

The lower panels of Figure 4 show the biases in measure- 
ments of the gas mass fractions. Since the hydrostatic mass 
is biased low, the derived / gas are typically biased high. The 
bias and scatter in / gas are especially large for the unrelaxed 
clusters, since the biases in M gas and M tot are added construc- 
tively. Both bias and scatter are significantly reduced for re- 
laxed clusters; for example, / gas (< ^sooc-* aie accurate t0 about 
10% at z — and 0.6. The biases become larger in cluster out- 
skirts, and the biases in / ga s(< r 5on c ) are about 18% at both low 
and high redshifts. 

In practice, additional biases in the estimated cluster masses 
(M est ) could arise from a bias in the estimation of a cluster 
virial radius. Figure 5 shows the biases in the estimated hy- 
drostatic mass within the estimated virial radius, M est (< r est ), 
relative to the true cluster mass, M true (< r true ), measured in 
simulations. The column indicated as M tot (< r est ) in Table 2 
summarizes average biases and scatter in the estimates of to- 
tal cluster masses computed this way. An underestimate of a 
cluster virial radius results in the increased bias in the M tot es- 
timate and an underestimate of the derived M gas by a similar 
amount, while leaving / gas relatively unchanged. It is these er- 
rors that contribute to the differences in the mass-temperature 
relations discussed in § 4.5. 

A similar study of biases in the hydrostatic M tot estimates 
has been done recently by Rasia et al. (2006). The cos- 
mological simulations were performed using the Gadget-2 
SPH code (Springel 2005) and the mock X-ray data were re- 
duced closely following the procedure of Ettori et al. (2002). 
Their sample included 5 clusters, including two with relaxed 
morphology, for which the comparison with our our results 
is most relevant. The mass biases in the most comparable 
case (direct hydrostatic estimates and reduced Chandra back- 
ground level) are -30% and -28% at ^25000 and -15% and 
-32% at r5ooc- These values are on the lower side of our dis- 
tribution for relaxed clusters (Fig.4). In particular, none of 
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Fig. 6. — Comparisons of various average ICM temperatures for the z = 
sample. The temperature averages are defined in the text. The filled and open 
symbols indicate clusters with relaxed and unrelaxed morphology, respec- 
tively. For each cluster results for the three projections are shown. 



our relaxed clusters shows such strong biases at r = r25ooc- 
However, it is difficult to compare our results directly. Not 
only the cosmological simulation codes are quite different, but 
also the data analysis algorithms are completely independent 
and significantly different. Investigations of the sources of 
discrepancy will require cross-checks (e.g., reduction of our 
mock data with the Rasia et al. pipeline) and/or using a larger 
sample of Gadget-2 clusters. 

4.4. Average Temperatures 

Average cluster temperature, (T), is another important ICM 
diagnostic and a key observable for cosmological application. 
Since the ICM is not isothermal within a cluster, the defi- 
nition of (T) is not unique. The differences among various 
definitions should be calibrated and taken into account. Fig- 
ure 6 compares different (T) obtained from the mock Chandra 
analysis, r spec , 7\pec r500c] > and T mg , as well as those measured 
directly from the 3D properties of gas in simulations, rj™ 
and T^ 1 , for all clusters at z — and 0.6 (see § 3.2). The gas 
with T < 0.086 keV (10 6 K) is excluded from all calculations 
of average temperatures, since it does not contribute to the X- 
ray flux. Note also that detectable small clumps are excluded 
from the mock Chandra analysis and hence do not affect the 
determination of (T). Average ratios of different temperature 
definitions are summarized in Table 3. 

The best agreement is between the true and X-ray derived 
gas mass-weighted temperatures. Other definitions show a 
constant offset relative to 7"™ with some scatter around the 



mean. For the relaxed clusters at z = 0, we find r£ M : T, 



: T, 



SIM 



spec 



1.19 : 1.13 : 0.99 : 1. The ratios are slightly 
different for the non-relaxed clusters, : r spe c : T m : 

= 0.99 : 1.07 : 0.96 : 1. The T spec is higher than T mg 
because the former is dominated by the inner, hotter region 
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TABLE 2 

Bias in M tot , M gas and / gas measurements. 











Bias ± Scatter 




Redshift 


Radial range 


Sample 




lvl wl ' est/ 


'"gas 


/gas 1 " 






all 


n 1 9 1 -t- o 1 ^6 

— U. 1^1 It U.I JO 


-0.222 ±0.181 


0.006 ± 0.047 


909 -t- 4^1 


2 = 


< >"2500c 


relaxed 


-0.075 ± 0.047 


-0.134 ±0.059 


-0.003 ± 0.017 


0.091 ± 0.062 






unrelaxed 


-0.146 ±0.163 


-0.270 ± 0.207 


0.010 ± 0.056 


0.268 ± 0.525 






all 


-0.163 ±0.095 


-0.253 ±0.162 


0.041 ± 0.048 


0.264 ±0.184 


2 = 


< ^(Mc 


relaxed 


-0.130 ±0.096 


-0.195 ±0.124 


0.011 ±0.023 


0.176 ±0.136 






unrelaxed 


-0.182 ±0.091 


-0.285 ±0.173 


0.058 ±0.051 


0.312 ±0.190 






all 


-0.122 ±0.116 


-0.245 ±0.151 


0.010 ± 0.045 


0.175 ±0.216 


2 = 0.6 


< r2500c 


relaxed 


-0.100 ±0.024 


-0.178 ±0.038 


0.000 ± 0.022 


0.112 ±0.038 






unrelaxed 


-0.135 ±0.145 


-0.286 ±0.178 


0.016 ± 0.054 


0.213 ± 0.266 






all 


-0.089 ± 0.206 


-0.144 ±0.306 


0.024 ± 0.072 


0.180 ±0.272 


2 = 0.6 


< rsooc 


relaxed 


-0.152 ±0.069 


-0.214 ± 0.079 


0.007 ± 0.029 


0.196 ±0.111 






unrelaxed 


-0.052 ± 0.250 


-0.102 ±0.379 


0.034 ± 0.087 


0.171 ±0.335 



"Estimated masses, M est , enclosed within r triie and r est are compared with true cluster masses, M trlle (< r tme ), 
measured in simulations. 
b Enclosed M gas within rsooc and r25ooc excluding the mass enclosed within 0.075r5ooc- 
C A fraction computed using above two quantities. 



TABLE 3 

Ratio of average cluster gas temperatures 









Average ± Scatter 




Redshift 


Sample 


7 1 /7^SIM 
1 mg/J mg 


7* /7^SIM 
1 spec/ 1 mg 


T SIM /T SIM 
f ew / 1 mg 


7^[0.5r500c] iji 
1 spec / 1 spec 


2=0 


all 

relaxed 
unrelaxed 


0.973 ± 0.040 
0.990 ± 0.029 
0.960 ± 0.042 


1.095 ±0.078 
1.126 ±0.049 
1.070 ±0.087 


1.078 ±0.283 
1.188 ±0.192 
0.992 ±0.315 


1.072 ±0.027 
1.070 ±0.015 

1.073 ±0.034 


2 = 0.6 


all 
relaxed 
unrelaxed 


0.985 ± 0.047 
0.985 ± 0.029 
0.985 ± 0.054 


1.080 ±0.075 
1.110 ± 0.056 
1.066 ±0.079 


1.511 ±0.619 
1.725 ± 0.789 
1.404 ±0.459 


1.044 ±0.019 
1.048 ±0.007 
1.042 ±0.022 



(see examples of the temperature profiles in Fig. 2 and 3). A 
tight correlation exits between r spec and r^ g M with the object- 
to-object scatter as small as 5% for the relaxed clusters and 
9% even for the unrelaxed systems. Although computed di- 
rectly from the simulation outputs, rj™ shows a much poorer 
correlation with T^ 1 than any of the observationally derived 
average temperatures. The scatter in the temperature ratios is 
generally higher for non-relaxed systems. 

The X-ray spectral temperature, r spec , is of particular inter- 
est to X-ray observers, because it is the most easily measured 
spectral characteristic. Our study shows that r spec is related 
to T mg by a constant factor, T sp<x /T mg « 1.14, for the re- 
laxed clusters and a slightly smaller factor, T spec /T mg = 1.12, 
in the unrelaxed systems. For each sub-sample, the ratio 
Tspec/T mg and its scatter are nearly identical at z = and 
0.6 and do not show any sign of evolution with redshift. Sim- 
ilar results are obtained for r s !p e ^ r50oJ , which is related to r spec 
by a constant factor, r t pe 5 c r50Oc] /T spe c ~ 1-04 - 1.07 with lit- 
tle redshift trend and scatter, indicating that r s ! pe ^ 500c] can be 
used as a reliable substitute for r spec . Finally, the ratio of 
T sp ec/T mg and its scatter from cosmological simulations is in 
good agreement with the Chandra results for low-redshift re- 
laxed clusters (Vikhlinin et al. 2006), where the average ratio 
in [0.15,l]r 50 oc is 1.15 and scatter is about 0.08. 



4.5. Implications for the Mass-Temperature relation 

The results in § 4.3-4.4 have direct implications for the 
mass-temperature (M - T) relation. The normalization of the 



M -T relation, M 5 , can be defined as 

E(z)M 5O0c = M 5 y^j (8) 

where E(z) = (Qm(1 + z) 3 + Qa) 1 ^ 2 for a flat universe with 
a cosmological constant. Table 4 lists the best-fit values of 
Ms and its uncertainties (arising because of the finite sample 
size) for different temperature averages (see also Table 2 and 
Fig. 6 in Kravtsov et al. 2006b). The values of Ms are com- 
puted using both the true mass, M l ^ c (< r™ ), measured in 
simulations and the estimated mass, M^^Xk r^ Qc ), derived 
from the mock analysis. For M^ c , the Ms is higher by about 
10% than their observed counterparts at both low and high- 
redshifts: M 5 = 3.06 + 0.16 (r spec ) and 3.64 ± 0.18 (T mg ) at 
z ~ (Vikhlinin et al. 2006) and M 5 = 3.36 + 0.32 (r spec ) 
at z ~ 0.5 (Kotov & Vikhlinin 2006). 7 The M 5 based on 
M™ , on the other hand, is lower by about 10%— 15%. In 
other words, the observed relations lie between the simula- 
tions derived relations based on the true and estimated masses. 

5. DISCUSSION AND CONCLUSIONS 

We present mock Chandra analyses of cosmological cluster 
simulations and assess X-ray measurements of galaxy cluster 
properties. To test observational X-ray procedures, we con- 
struct mock Chandra images of the simulated clusters and 

7 To compare the M-T relations from simulations and observations in 
the same radial range (0. 15f5ooc-''500c). we applied the correction factors of 
(T(0. 15r 50 oc - r 5Q0c )/W0 kpc - r 500 c)> = 0.97 (for 7/ spec ) and 0.94 (for 7/ mg ) 
to the observed relations. 
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TABLE 4 

Normalizations of the Msqq c - T relation 









m s , io 14 r ] M 


Mass" 


Redshifi 


Sample 




^spec 


j,[0.5r5ooc] 
1 spec 


Mtrue 
m 500c 


2 = 


all 
relaxed 
unrelaxed 


4.46 ±0.18 
4.02 ±0.10 
4.80 ±0.19 


3.78 ± 0.21 
3.32 ±0.10 
4.14 ±0.22 


3.40 ±0.17 
3.00 ±0.09 
3.71 ±0.18 


M ,ruc 


z = 0.6 


all 

relaxed 
unrelaxed 


4.65 ±0.17 
4.50 ±0.12 
4.73 ±0.18 


4.07 ±0.17 
3.77 ±0.12 
4.22 ±0.17 


3.81 ±0.16 
3.51 ±0.12 
3.96 ±0.16 


M 500c 


2 = 


relaxed 


3.14 ±0.03 


2.59 ± 0.03 


2.35 ± 0.04 


M est 
m S00c 


2 = 0.6 


relaxed 


3.55 ± 0.07 


2.97 ± 0.06 


2.77 ± 0.06 



a tf 5 ™ is the true M tot of a cluster enclosed within the true r^, measured di- 
rectly in simulations; M^ 0c is the estimated hydrostatic mass enclosed within 
the estimated rSL from the mock data analysis. 



derive properties of X-ray clusters using a model and proce- 
dure essentially identical to those used in real data analysis. 
The sample includes 16 clusters spanning a broad mass range 
(5 x 10 13 - 2 x 10 15 A -1 M ) simulated in the ACDM cosmol- 
ogy with high spatial resolution and including various physi- 
cal processes of galaxy formation, such as radiative cooling, 
star formation and other processes accompanying galaxy for- 
mation. We analyze the simulated clusters at z = and their 
most massive progenitors at z = 0.6. To test the sensitivity 
of the X-ray analysis to the cluster dynamical state and sub- 
structure, we also distinguish unrelaxed clusters from relaxed 
systems based on the overall structural morphology of the 
Chandra images, as usually done to classify observed clus- 
ters. Below we discuss and summarize our main findings and 
conclusions. 

5.1. Reconstruction of 3D ICM Properties 

Derivation of the cluster 3D properties from projected X- 
ray data usually relies on simplifying assumptions such as 
spherical symmetry of a cluster or small-sale uniformity of 
the ICM. Additional complications arise from a non-ideal re- 
sponse of X-ray telescopes and limited statistical quality of 
data. This leads to several effects which are hard to reproduce 
analytically without creating mock X-ray data from realistic 
cluster simulations. For example, clumps of colder gas as- 
sociated with massive subhalos can be detected and masked 
out from a further analysis while smaller subhalos remain un- 
detected and bias projected temperature and X-ray brightness 
measurements (Mazzotta et al. 2004; Vikhlinin 2006). One 
of the goals of our work is to include all such effects as com- 
pletely as possible in generation of the mock X-ray data and 
thus to test how they bias X-ray measurements and recovery 
of 3D ICM profiles. 

Reassuringly, we find that the ICM models and analysis 
method presented in Vikhlinin et al. (2006) provide a good de- 
scription of emission measure and projected temperature pro- 
files of the ICM. For relaxed clusters, these models recover 
the 3D gas density, temperature, and gas mass profiles with 
the accuracy of a few percent at 0.1 < r/>5oo c < 2.0. The 
models are flexible enough to describe also non-relaxed sys- 
tems and provide reasonably accurate gas mass and temper- 
ature measurements. A practical implication of these results 
is that projection effects, cluster substructure, and deviations 
from spherical symmetry do not strongly affect reconstruc- 
tion of 3D ICM density, temperature, pressure, and entropy 
profiles from X-ray data. Therefore, observational results on 



ICM profiles can be directly compared to the properties of 
simulated clusters without a need for detailed mock X-ray 
analyses. Such a comparison for our sample will be presented 
in the future paper (Nagai et al. 2006). 

The total ICM mass is measured quite accurately in all our 
simulated clusters. Results are most accurate for relaxed sys- 
tems: biases and scatter in M gas measurements are accurate to 
better than ~ 3% percent, independent of a cluster redshift or 
a radial range within in which the M gas is measured. For un- 
relaxed systems, the scatter is larger, as expected (Mathiesen 
et al. 1999), but the bias is still small (< 10%). 

Average temperature is an important diagnostic of proper- 
ties and physical processes of the ICM and a key observable 
for cosmological applications. Since clusters are not isother- 
mal, the definition of the (T) is not unique nor can (T) be 
completely accurately derived from data. Mock X-ray anal- 
ysis is therefore required to study the relation among differ- 
ent (T) definitions used in observational and theoretical stud- 
ies. Our analysis is focused on the comparisons of the true 
gas-mass-weighted average temperature, T^ 1 , and the X-ray 
spectroscopic r spec or gas-mass-weighted T mg temperatures 
derived from the mock data. We find that these temperatures 
are not identical, but the difference is mostly a constant factor 
with a relatively small object-to-object scatter. For example, 
T S pec/T mg ~ 1.14 ± 0.05 for relaxed clusters, while for unre- 
laxed objects this ratio is T spec /T mg « 1.12 + 0.09. All the 
(T) definitions we tested are therefore equivalent, but the dif- 
ference should be kept in mind when, e.g., observed M - T 
relations are compared with results of simulations. A notable 
exception is the "emission-weighted" temperature defined as 
r ew = fTp 2 T l/2 dV/fp 2 T 1/2 dV (Bryan & Norman 1998; 
Frenk et al. 1999; Muanwong et al. 2001). Even computed di- 
rectly from the simulation outputs, this quantity shows a much 
poorer correlation with 7"™ than any of the observationally 
derived average temperatures (Fig. 6). Our results thus sup- 
port the suggestion of Mazzotta et al. (2004) that T ew should 
not be used as a measure of the average cluster temperature. 

5.2. Accuracy of Hydrostatic Mass Estimates 

One of the major applications of the gas density and temper- 
ature profile measurements from the X-ray data is the estimate 
of a gravitationally bound mass of a cluster assuming that the 
ICM is in hydrostatic equilibrium in the cluster potential. The 
X-ray method is by far the most precise way to estimate M tot 
in individual objects at large radii. This limits our ability to 
test biases in the X-ray derived M tot estimates through com- 
parisons with independent observational techniques such as 
weak lensing. The analysis of our mock X-ray data, however, 
provides important clues. 

We find that the hydrostatic mass estimates from our mock 
data are biased low by about 5%-20% throughout the virial 
region (Fig.4), even in clusters that are identified as relaxed 
by their X-ray morphology. Such biases seem to be too high 
to be attributed to inaccuracies in 3D gas density and temper- 
ature profile reconstruction. A more plausible explanation is 
a departure of the ICM from a hydrostatic state. Indeed, an 
additional pressure component due to residual random bulk 
motions of the ICM in our simulated clusters completely ex- 
plains the bias in the hydrostatic mass estimates (this topic 
will be discussed in greater detail in Lau et al. 2006). Sim- 
ilar results have been also obtained by independent simula- 
tions performed with different codes and implementation of 
physical processes (Rasia et al. 2006; Dolag et al. 2005). The 
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M lot bias derived in the mock X-ray analysis of our relaxed 
clusters increases towards their outskirts. The average levels 
are -8% at r = r2500c and -13% at r = rsooc, and quickly 
become very large beyond r^ooc (Fig. 2), which qualitatively 
mimic the radial dependence of turbulent pressure observed 
in previous simulations (e.g., Evrard et al. 1996). In unre- 
laxed clusters, not only the departures from hydrostatic equi- 
librium are strong but also the density and temperature gra- 
dients are not measured accurately (Fig. 3) resulting in M tot 
biases of 20% or more at all radii. 

It is unclear to what degree the level of ICM turbulence 
found in these simulations applies to real clusters because 
most of the codes model ICM as an ideal inviscid fluid with 
some amount of numerical viscosity. Further progress in this 
area using cosmological simulations will have to rely on bet- 
ter understanding and modeling of viscosity of the weakly 
magnetized plasma (Ruszkowski et al. 2004a,b; Dolag et al. 
2005; Sijacki & Springel 2006). Direct measurements of 
the ICM turbulence through Doppler broadening of emission 
lines (Inogamov & Sunyaev 2003) must await a launch of X- 
ray calorimeters, although some clues could be provided also 
by observations of X-ray brightness fluctuations (Schuecker 
et al. 2004). Keeping these uncertainties in mind, we can treat 
the M to t biases quoted above as upper limits. 

Cosmic rays and magnetic fields may also provide non- 
thermal pressure support which is not accounted for in the 
hydrostatic estimates nor included in our simulations. In addi- 
tion, temperatures of electron and ion components of the ICM 
are not necessarily the same, because strong shocks heat pri- 
marily ions and the collisional electron-ion equilibration time 
is long in the cluster outskirts (Markevitch et al. 1996; Fox 
& Loeb 1997). If T e (derived from X-ray spectra) is indeed 
lower that T t (unobserved), the hydrostatic mass estimates are 
biased low. We note, however, that at least in one case there is 
evidence for fast equilibration of the electron and ion temper- 
atures in a post-shock region of the ICM (Markevitch 2005). 

While direct measurements or theoretical modeling of the 
non-thermal pressure support may be difficult, it can be con- 
strained by comparisons of the hydrostatic mass estimates 
with those derived from gravitational lensing for the same ob- 
jects. An important caveat is that lensing measures the mass 
in the projected aperture while the mass within the sphere 
of the same radius is rather uncertain because of projection 
effects. Metzler et al. (2001) estimate that Msm c cannot be 
determined from weak lensing with better than ±30% uncer- 
tainty. Any comparisons of X-ray hydrostatic and weak lens- 
ing masses must be done in an average sense, e.g. through de- 
rived normalizations of the correlation between the M tot and a 
low-scatter mass proxy such as Y x (Kravtsov et al. 2006b) or 

^spec- 

To summarize, our results indicate that an X-ray analysis 
of relaxed clusters correctly recovers gas density and thermal 
electron pressure gradients. The simulations provide an upper 
limit (5%-20% depending on radius) for the mass bias due to 
turbulent motions. The contribution of other pressure compo- 
nents will have to be constrained by future observations. 

5.3. Implications for Cluster-Based Cosmological Tests 
5.3.1. Tests Based on / gas 

X-ray measurements of the baryonic mass fraction in mas- 
sive clusters can be used to constrain QjfOr/i (White et al. 
1993) and also provide a potential standard candle (Sasaki 
1996; Pen 1997; Allen et al. 2004). The major observational 



ingredient in these tests is a measurement of the hot gas mass 
fraction within a sufficiently large radius. Our analysis pro- 
vides clues for biases in / gas measurements from high-quality 
Chandra data. We find that the / gas determinations are bi- 
ased high, primarily because the hydrostatic method underes- 
timates the total cluster mass. As was discussed above, the 
bias in the hydrostatic estimates is related to the physical pro- 
cesses in the ICM (e.g., turbulence, cosmic ray and magnetic 
pressures etc.) and not to the deficiencies of the X-ray anal- 
ysis. Further progress in confirming validity of / gas measure- 
ments will have to rely on independent determinations of clus- 
ter total mass (see discussion in § 5.2). Some sources of the 
bias (e.g., turbulence) could be minimized if we focus only 
on relaxed systems. We find that in such systems, both bias 
and scatter in / gas can be controlled to within about 5%-8% at 
72500c and 15%-20% at r 50 o c , which can be considered the up- 
per limits on the / gas measurement bias caused by turbulence 
(see § 5.2). 

Note that if the ICM turbulence and other processes dis- 
cussed in § 5.2 do play a role in biasing the hydrostatic mass 
estimates, the bias can potentially be redshift-dependent since 
the merger rate increases at higher z (e.g. Gottlober et al. 
2001). We do not find any detectable redshift dependence of 
the bias in the / gas measurements in relaxed systems. How- 
ever, our upper limits, ~ 5%, are weaker than the accuracy 
required, e.g., for using / gas as a standard ruler in the dark en- 
ergy studies (Allen et al. 2004). Mock analysis of much larger 
cluster samples will be required to properly address this issue. 

5.3.2. Tests Based on the Cluster Mass Function 

Another group of cosmological tests are those based on 
measurements of the cluster mass function. The relevant is- 
sue for our work is how well various mass proxies can be 
measured. A related question is how well the cluster mass vs. 
proxy relation can be calibrated by the X-ray analysis, which 
was discussed in § 5.2. 

We did not discuss the easiest X-ray mass proxy, the to- 
tal luminosity, because it is also the least accurate and least 
reliably reproduced by numerical simulations. The next, in 
terms of being easily measured, is the average ICM tempera- 
ture. The practical averages to use are r spec or TspeJO.Srsooc] 8 - 
These definitions are freely interchangeable since the ratio 
rspec[0.5r5ooc]/r s pec is nearly the same in all clusters (Fig. 6, 
Table 3). One disadvantage of the spectral temperature as 
a mass proxy is that its measurement is somewhat sensitive 
to the cluster dynamical state and small-scale substructure. 
This sensitivity arises because r spe c is rather sensitive to the 
presence of lower-temperature components. This effect con- 
tributes to the systematic difference in the normalizations of 
the M - T relation between relaxed and non-relaxed clus- 
ters and also slight deviations from the self-similar evolution 
(Table 4). Some implications of this effect for the cosmo- 
logical measurements are discussed by Rasia et al. (2005). 
The prospects for calibrating the r spec measurement biases by 
cosmological numerical simulations are mixed. The effects 
related to the cluster dynamics, such as merger activity and 
shocks etc., are treated sufficiently accurately. However, the 
biases caused by the presence of undetected cold clumps will 
require a more realistic modeling of the ICM cooling and as- 
sociated feedback. 

More accurate X-ray proxies for cluster total mass is the 

8 Other definitions such as r mg require measurements of the temperature 
profiles which is feasible only for a few best-observed clusters. 
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mass of the hot gas, M gas , and the "X-ray F-parameter", 
Yx = M gas x r sp ec (Kravtsov et al. 2006b). The properties of 
Y x (this is the best X-ray mass proxy known) are extensively 
discussed elsewhere (Kravtsov et al. 2006b), so we concen- 
trate on M gas here. The M gas is measured very reliably. The 
average bias is < 5% with only a few strong outliers (Fig. 4) 
and there is no detectable redshift dependence. A similar ac- 
curacy should in principle be achieved for the M tot estimates 
since to the first order, M tot cc M gns . Note that the biases 
in Mgas measurements are mostly due to large-scale devia- 
tions of the cluster body from spherical symmetry. This effect 
should be reliably treated in the current generation cosmolog- 
ical simulations, and so the accuracy of the M g a S -based total 
mass estimates can be improved yet further. All we need is a 
larger sample of simulated clusters. A major obstacle, how- 
ever, is the uncertainties in the mass fraction of cold baryons 
(stars and molecular gas). This fraction can be substantial, as 
indicated by our simulations (Kravtsov et al. 2005) and sig- 
nificant deviations of the observed / gas in clusters from the 
cosmic mean (e.g., Ettori 2003; Vikhlinin et al. 2006). Fur- 
ther progress in this area requires significant improvements 
in the reliability of the stellar mass measurements and treat- 
ment of the star formation in the simulations. We will further 
discuss these issues in the forthcoming papers. 

To briefly summarize our results, we find that the analy- 
sis of deep X-ray observations of galaxy clusters can reliably 
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recover its intended targets: the distribution of density and 
temperature of the hot ICM. The accuracy of estimating the 
total cluster mass from these parameters is limited by addi- 
tional physical processes in the ICM, such as turbulence, mag- 
netic and cosmic ray pressure, and possible departures from 
the electron-ion equilibrium. Our work provides a realistic 
estimate of the effects of turbulence, which are found to cause 
underestimation of M tot by 5%-20% in clusters that can be vi- 
sually classified as relaxed (this can be considered as an upper 
limit for the bias if the ICM viscosity is non-negligible). Con- 
straining the role of other effects will probably have to rely 
on accurate independent mass measurements (e.g., by weak 
lensing) in a large, complete cluster sample. 
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